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j3 ■ Abstract 

^ I We study a recent conjecture about the behavior of the quantum relative entropy 

compared to the relative entropy of entanglement in open bipartite systems. The con- 
jecture states that, under a dissipative time-evolution, the positive rate of change of the 
relative entropy will always be larger than that of the relative entropy of entanglement. 



' After explicitly solving a two-qubit master equation of Lindblad-type with separable 

and entangled stationary states, 
states with an entangled asympt 
asymptotic states are separable. 



On ■ and entangled stationary states, we show that the conjecture can be violated for initial 

^ states with an entangled asymptotic state, while it appears to be confirmed when the 



O 
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: 1 Introduction 

> ■ 

^ ■ The importance of quantum entanglement as a physical resource for performing informa- 

■ tional tasks which would be classically impossible [1] has spurred the study of its dynamical 
behavior in many different systems. The time-evolution of most of these is reversible and 
generated by a Hamiltonian; however, it is important for quantum entanglement to be used 
as an efficient physical resource that its temporal behavior should also be studied when 
systems are driven by noisy environments and their dynamics is irreversible. 

In the following, we will consider open quantum systems [21 |3l H] i.e. systems where 
the interactions between the subsystem S and the external environment E, though weak, 
cannot be neglected. In this case, a standard way of obtaining a manageable dissipative 
time-evolution of the density matrix gt describing the state of S at time t is to construct 
it as the solution of a Liouville-type master equation dtQt = L[f?t], where the generator L 
of Lindblad type [5l |6] takes care of the effects of the environment through a characteristic 
matrix of coefficients known as Kossakowski matrix. This can be done by tracing away 
the environment degrees of freedom and by performing a Markovian approximation, i.e. by 
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studying the evolution on a slow time-scale and neglecting fast decaying memory effects. 
Then, the irreversible reduced dynamics of S is described by one-parameter semigroups of 
linear maps, called quantum dynamical semigroups, obtained by exponentiation: 7^ = e*^, 
t > 0, such that gt = 7t[f?]- In order to guarantee full physical consistency, namely that 
id ® 7t be positivity preserving on all states of the compound system S + Sd for any inert 
ancilla Sd, 7* must be completely positive [2l[3|,[7]. 

The formalism of open quantum systems has been used to describe the tendency to thermal 
equilibrium of a small system in weak interaction with a large heat bath at a certain temper- 
ature. The main tool in this thermodynamical picture is the quantum relative entropy [3]; it 
is related to the difference between the free energy of the irreversibly evolving open quantum 
system and that of its equilibrium asymptotic state: this difference monotonically decreases 
in time because so does the quantum relative entropy with respect to completely positive 
maps |8] , as quantum dynamical semigroups are. Namely, the time derivative of the quantum 
relative entropy, called entropy rate, has a definite sign. 

The quantum relative entropy has also been used as a possible measure of the entanglement 
content of a quantum state: the so-called relative entropy of entanglement provides a pseudo- 
distance between a state and the closed convex set of separable states [9]. 

In pD] the natural question was raised whether the entropy production due to thermo- 
dynamical tendency to equilibrium is somewhat related to the entanglement rate, that is to 
the speed of variation of the relative entropy of entanglement. A conjecture was put forward 
that for systems immersed in an external bath without a direct source of entanglement due 
to Hamiltonian interactions, the absolute value of the entanglement production is always 
smaller than the entropy production. 

Typically, a system S immersed in a large environment E is subjected to decoherence; 
therefore, one expects quantum entanglement to be generically depleted by a dissipative and 
noisy time-evolution. The conjecture mentioned above is motivated by the fact that, if a 
quantum open system tends to a separable equilibrium state, then, in a suitable neighborhood 
of the latter, the entanglement production is zero while the entropy production is not. Indeed, 
in [To] a concrete example that validates the conjecture is offered of a 2-qubit system in which 
only one of them evolves as a quantum open system. In such a case an initial maximally 
entangled state evolves towards a separable steady state with an entropy production always 
larger than the speed with which entanglement is dissipated. 

However, in certain specific situations, an environment affecting both parties of a bipar- 
tite system may even build quantum correlations between the subsystems which compose 
S (see, e.g., [HI [121 1131 13 [II])- Iii particular, in [7] this possibility is shown to depend on 
the specific form of the generator of the reduced dynamics. In [TS] an inequality was found, 
involving the entries of such a matrix which, if fulfilled, is sufficient to ensure that a specific 
initial separable pure state of two qubits gets entangled. Further, in [16] this inequality was 
proven to be a necessary and sufficient condition for environment-induced entanglement in 
an initially separable pure state of two qubits. Even more interestingly, starting from an ini- 
tially separable state, the entanglement generated at small times can persist asymptotically; 
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also, starting from an initially entangled state, its entanglement content can asymptotically 
increase. 

This work is organized as follows: in Section 2, we consider the open dynamics of two 
qubits with a generator that depends on a parameter which allows to range over all the above 
mentioned cases, analytically solving the master equation; then, in Section 3, we introduce 
the notions of entropy and entanglement rates and the conjecture from [10]; finally, in Section 
4, numerically studying the time-behaviour of the entanglement and entropy rates for various 
initial states, we show that, whenever there is asymptotic entanglement the conjecture in \XOl 
is violated, while it holds if there is no asymptotic entanglement. 



2 The Reduced Dynamics 

Let a bipartite system composed of two qubits be immersed in an external environment in 
such a way that, via standard weak-coupling limit techniques [2], one describes their reduced, 
irreversible dynamics by means of the master equation 

dtQt = L[gt] = -i — Qt + (^Sj^Sj - -{S^S^, , (1) 

where fl is the system frequency, := cTj (g) I + I (g> cTj, I is the 2x2 identity matrix, a^, 
i = 1,2,3 are the Pauli matrices and the matrix 




A = [Aij] = 1 , a e M , < 1 , (2) 



is positive semi-definite. This latter request ensures that the semigroup generated by ([I]) 
consist of completely positive maps '-ft for alH > [2]. 

Remark 1 By means of the single qubit Pauli matrices a^^^ = cTj (8> I and a^'^^ = I (g) (jj one 
writes the purely dissipative contribution to the generator as [7j 

DIft] = f: A, f {"l" p. "f - \ {-f'-.'"' , p.}) . (3) 

i,j=l a,b=l 

In this way there are six Kraus operators (tI"'\ a = 1,2, i = 1,2,3 and the 6x6 Kossakowski 
matrix reads 

From the theory of open quantum systems [21 [5l [3l H] one knows that the coefficients K^j^^ 
in the Kossakowski matrix relative to the i-th Pauli matrix of the a-th qubit, respectively 
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the j-th Pauli matrix of the 6-th qubit, a,b = 1,2, i,j = 1,2,3, are determined by the 
Fourier transforms of the two-point time-correlation functions with respect to an environment 
equihbrium state oo, uj{bI_"'^ B'j'\t)), of the environment operators B^""^ appearing in the 

system-environment interaction Hj = Yli=i(^'^i^'' ® -B^^"* + o"!^-* (g) b'^'^^^. The symmetric 
form of (|3]) thus results when both qubits are linearly coupled to bath operators such that: 
^S,3 = ^Sa = ^1,2,3 and uiB^^^B^it)) = 0. 

Remark 2 Considering two qubits weakly interacting with a thermal bath modeled as a 
collection of spinless, massless scalar fields (see, e.g., [7]) at very high temperature T = 1//3, 
the parameter a in the Kossakowski matrix is related to P, i.e. a = —f3^l, where Q is the 
system frequency when isolated from the environment. Correspondingly, in the case of one 
qubit immersed in such a thermal bath at high temperature (/3 ^ 1), any initial state is 
driven to the thermal asymptotic state 

exp(-/3fia3) 1/ 
2cosh/31] "2 1^-^^"^ 



The master equation ([T]) is explicitly integrated in Appendix A; in the following we will 
mainly focus upon the time-evolution of initial states of the form 

^ = a|l)(l| +rf|2)(2| + 6|3)(3| +c|4)(4| , a,b, c, d e , a + b + c + d = 1 , (5) 

diagonal with respect to the orthonormal vectors 



|1) = |00), |2) = |11), |3) 



|01) + |10) ^ |01)-|10) 



V2 



V2 



(6) 



where cxslO) = |0), o"3|l) = — 11), and |00), |01), |10), |11) form the so-called standard basis 
in ® C^, with respect to which the states ([5]) are represented by 



fa 0\ 

^ ^ 

^ ^ 

\0 dj 



(7) 



From equations flM|) -f p9|) and fHT]) in Appendix A, it turns out that these initial states evolve 
at time t > into states of the same form 



gt = at\l){l\ + dt\2){2\ + bt\3){3\ + q|4)(4| , 



(8) 
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where Cf = c and 



at 



dt 



+ 



+ 



R + VT^ 



{l + afa -2{l-a)d + {l + afb 



3 + a2 

2(1 + a) a - (1 + 



(l + a)(3 + a2) 



3 + a' 



R - VT^ 



a- 



2(1 + a) a - {l-af{b + d) 



3 + a2 (l-a)(3 + a2) 

(l + a)2a -2{l + a)d +{l + ayb 



E-it) 



3 + a2 



E+{t) 



(l-a2) 



R + yr^ 



(1 + afa+{l- af - 2(1 - a^) 5 



3 + a2 (3 + a2)(i_^2) 

2(1+^2)^ _ (1 _^2)(^^^) 



E-{t) 



3 + a2 

with i? = a + 6 + rf = l — c and 



E+{t) 



(9) 



(10) 



(11) 



^_^(t) = e"^* cosh4tVl - «2 , E^{t) = e"^* sinh4tvr 



Since hmf_>.+oo -E'±(^) = 0, the asymptotic states resulting from the initial states ([5]) are 
Qooic) = (1 - c) |1)(1| + (1 - c) |2)(2| + (1 - c) |3)(3| + c |4)(4| . 



3 + a2 



3 + a2 



3 + a2 



(12) 

There is thus a one-parameter family {goo{c)}o<c<i of asymptotic states such that all initial 
states of the form ([5]) with the same c go into the same Qooic). 

In order to study the asymptotic entanglement generation capability of the present model, 
we shall measure the entanglement of 2-qubit states p by the concurrence |19j : 

C{g) = max{0, Ai - A2 - A3 - A4} , 

where Ai > A2 > A3 > A4 > are the square roots of the positive eigenvalues of gg with 
^ = o"2 ® f?* (72 ® 0-2, g* denoting the complex conjugated matrix. For 2-qubit states of the 
form ([7]), p = p* and one easily computes 



/do 0\ 

^ ^ 

^ ^ 

\0 a/ 



gg 



fad 


Vo 



2 







adj 



This latter matrix has positive eigenvalues ad (twice degenerate), ^2^ c2; then, their square 
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roots Ai > A2 > A3 > A4 in decreasing order yield 
C{g) 



max 0, 2 ( : — ^ — Vad 



C(f?oo) 



max < 0, 



2 

1 - - 4c 



3 + a2 



(13) 
(14) 



Remark 3 As already emphasized in the introduction, despite decoherence, the presence 
of an environment need not have only destructive effects in relation to entanglement: entan- 
glement can even be asymptotically increased with respect to the initial amount. This can 
happen in the present case and the entanglement generation capability of the environment is 
entirely due to the non-Hamiltonian contribution ([3]) to the generator in ([1]). Indeed, the 2- 
qubit Hamiltonian does not contain coupling terms and cannot be a source of entanglement; 
instead, this can be true for because the off-diagonal contributions in the Kossakowski 
matrix (jl]) couple the two qubits. Of course, this is only necessary, but not sufficient to 
ensure entanglement generation and its asymptotic persistence. They indeed depend on a 
trade-off between the off-diagonal couplings and the purely decohering diagonal terms in 



3 Entropy and Entanglement Rates 

In this section we shall introduce the notions of entropy and entanglement rates; for sake of 
simplicity, we shall consider finite d-\evel systems whose states are described by normalized, 
positive, dx d density matrices g G M^(C). Given two such density matrices, their quantum 
relative entropy is defined by [8] 

5'(pi||p2) = Tr(^^i(log^i - log^2)) • (15) 

Consider a quantum system with Hamiltonian H; if in contact with a heat bath at tempera- 
ture T = 1//3 (with the Boltzmann constant k = 1), it is expected to be driven asymptotically 
into the thermal (Gibbs) state gx = e~^^ /Zp, where = Tr[e~^-'^]. Suppose that under an 
irreversible time-evolution g gt, a.n initial state g is driven into thermal equilibrium, that 
is limf^+oo gt = Qt\ then, 

^^(^?t||^T) = ^TT{gt{\oggt + \ogZp + PH)) = -T S{gt) + TrigtH) +T \ogZp , 

where S{g) = — Tr^plogf) is the von Neumann entropy of the state p. Since the second term 
corresponds to the system's internal energy, the first two contributions give the system's free 
energy corresponding to the time-evolving state pt [3]: 

F{gt) = U{gt) - T S{gt) , U{gt) = Tr(ft H) . 
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Finally, F^qt) = — \ogZp implies that the quantum relative entropy is related to the differ- 
ence of free energies 

S{Qt\\QT) = (^(F{Qt) - F{qt)) . 

Because of the second law of thermodynamics, the above quantity should be positive and 
its time-derivative non-positive. The first property is guaranteed by the properties of the 
quantum relative entropy |8j, while the second one holds true when the irreversible time- 
evolution is given by a Markovian semigroup, that is when Qt = ■jtiol and 7407^ = 7^07^ = ■j^.^^. 
for all s,t>0. Indeed, since ■jtlQx] = Qt, one derives 

S{gt\\gT) = s(^-ft[Q]\\it[QT]^ = s(^-ft^sOis[0]\\it~sOis[QT]^ 

< Si^s[Q]\\ls[QT]) = Sigs\\gT) VO<s<t, 

where the last inequality follows from the fact that the quantum relative entropy decreases 
under the action of completely positive trace-preserving maps |8]. 

Based on the previous thermodynamical arguments, one may consider generic open quan- 
tum dynamics g 1— )■ 7t[^] = Qt with asymptotic states limt_>+oo Qt = Qoo, that are not nec- 
essarily thermal ones. The speed of convergence to such stationary states starting from an 
initial state g will then be measured by the entropy rate 

cr[0t] = -^S{gt\\goo) = Tr (^^^(log ^00 - log ft)) . (16) 

The entropy production that accompanies the tendency to equilibrium of the states of the 
form (IH]) is easily computed; indeed, being the states gt and ^00 diagonal with respect to the 
same orthonormal basis, the entropy rate ( fT6i) has the analytic expression 

(l_«)2(i_c) . , (i-a^)(i-c) . , (l + a)2(i_c) 
aj(3 + a^) bt{3 + a'') dt[3 + a'^) 

When the density matrix p is the state of, say, a bipartite quantum system, it makes sense 
to introduce the relative entropy of entanglement, 

E[g] = inf S{g\\gsep) , (18) 

Qsep 

as a measure of the entanglement content of p. Indeed, the above quantity vanishes if and 
only if p is separable and can be used to measure the distance of p from the convex set 
of separable states; furthermore, it cannot increase, but at most remain constant, under the 
action of local operations, described by trace-preserving completely positive maps acting 
independently on the two parties [171 flS] . 

Analogously to what was done for the entropy production, one may look at the entangle- 
ment rate when the system evolves, i.e. at the time-derivative of the pseudo-distance 

aE[gt] = j^E[gt] . (19) 



^The relative entropy of entanglement is not exactly a distance since it is not symmetric. 
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In [To] it was argued that 

(rE[Qt]\<(y[Qt] (20) 

always holds in absence of direct entangling interactions between the parties. The argument 
on which the conjecture is based is that decoherence is expected to deplete entanglement 
before reaching the asymptotic state and thus before the entropy production vanishes. Such 
asymptotic intuition is then extrapolated at all times. 

Now we will illustrate the various possibilities offered by the reduced dynamics discussed 
in the previous section; in particular, we will compare the entropy and entanglement rates, 
thus checking the validity of the conjecture (I2U]) . 

We will first derive an explicit expression for the relative entropy of entanglement f|T8|) 
in the case of states of the form (|8]) and then compute numerically the behaviour of its 
time-derivative ffTI?]) . 

Let us first rewrite (fT8|) as follows: 

E{Qt) = -S{gt) - supTilgtlog Qsep) , (21) 

where S{gt) is the von Neumann entropy of the time-evolving state. The following result, 
which we prove in Appendix B, helps to explicitly solve the above maximization problem. 

Proposition In the case of states as in ([S]), the supremum in (1211) is achieved for separable 
states of the form 

gsep = x\l){l\ + u\3){3\ + v\A){A\ + y\2){2\ , (22) 
where the parameters x, y, u, v are real and such that 

\u — v\ 

x + u + v + y = l, — - — < ^/xy . (23) 
This leads to the following maximization problem 

Eigt) = -S{gt) - sup (atlogx + dthgy + bthgu + Ctlogvj , (24) 

which can be analytically solved. 



4 Results 

The above maximization problem is explicitly solved in Appendix B thus permitting to 
calculate numerically the entanglement rate (HM . to compare it with the entropy rate flTBl) 
and to check the conjecture (!20|) . We shall do this in a number of cases that cover all possible 
initial and asymptotic entanglement properties for which we plot the behaviors of the relative 
entropy and of the relative entropy of entanglement, separately, while the entropy rate and 
the entanglement rate are plotted together for direct comparison. 
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In the following, the choice of the range of values for the plots' axes was made only for 
graphic reasons to make the plots clearer. Moreover, in the first four cases, we took the 
parameter in the matrix A to he a = 0.5 as it makes the plots easier to read; changing a 
does not alter the results. In the last example, instead, we show two different behaviors of 
the entanglement of the initial state depending on the choice of the parameter a. 
Case 1. An initial pure separable state ( l25l) goes into a mixed separable state; the dis- 
sipative time-evolution is not able to generate entanglement at any time, as shown by the 
second and the third plot below where the entropy of entanglement and the entanglement 
rate are both zero. In this case the conjecture (!20!) holds. 

0=\1){1\ (25) 




0.1 0.2 0.3 0.4 0.5 -1.0 - 



Figure 1: Case 1: a = 0.5; Left: S{gt)] Right: E{gt) 

10 - \^ 
8 - \ 
6 - \ 
4 - 

2 - ""-^ 



Figure 2: Case 1: a = 0.5; cr[gt] dashed line, |o'£;[f)t]| continuous line 
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Case 2. An initial mixed separable state (126|) goes into a mixed entangled state and the 
conjecture (12UI) is violated after some time. 



^=i|3)(3| + ^|4)(4| 



(26) 



2(3 + a2 
C(f5)=0, C{q, 



(l + g)^ 
2(3 + a2 
2^2 



3 + a2 



> Va G [-1, 1] . 




Figure 3: Case 2: a = 0.5; Left: S{Qt\\Qoo)] Right: E[Qt 



0.30 
0.25 
0.20 
0.15 
0.10 
0.05 




Figure 4: Case 2: a = 0.5; a[Qt\ dashed line, |cT£'[^t]| continuous line 

Case 3. An initial mixed entangled state (!27j) goes into an asymptotic mixed state which 
is more or equally entangled and the conjecture (!20|) is violated after some time. 



^ = ^|1)(1| + — |2)(2| + — |3)(3| + — |4)(4| 
^ 10 ' ' 10 ' ' 10 ' ' 10 ' ' 



(27) 



2 3 + 4a=^ 2 
5 3 + a2 - 5 



Va G [-1, 1] 
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Figure 5: Case 3: a = 0.5; Left: S{gt\\Qoo)', Right: E[gt] 




0.2 0.4 0.6 0.8 



Figure 6: Case 3: a = 0.5; <j[gt] dashed hue, |cr£;[f)t]| continuous 



hne 



Case 4. An initial mixed entangled state ( l28l) goes into a state with less entanglement 

(2^ 



^=^|2)(2| + ^|3)(3| + ^|4)(4| 



3(1 -a)\ , , , 



5(3 + 



3(1 + af , , 3(l-a2) , , 2 ,,,,,, 

C(f?oo) = fora2<A 
oo. = Sfe^ < I for ^ < < 1 • 



With the choice a = 0.5, the dissipative time-evolution shows a sudden death of entan- 
glement, that is the concurrenc^ (fT3i) vanishes at finite time. The conjecture ( l20l) always 
holds. 



^In Figure 7, instead of plotting the concurrence as defined in (fT3|) . we have plotted the difference 
|6 — c| — 2Vad: this simply means that, as soon as this difference becomes negative, the state is separable. 
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0.30 - 
0.25 |- 
0.20 - 
0.15 - 



Figure 7: Case 4: a = 0.5; C{gt) 




Figure 8: Case 4: a = 0.5; Left: S{gt\\goo)] Right: E[gt 



Figure 9: Case 4: a = 0.5; cT[gt] dashed hue, continuous hne 



Case 5. An initial mixed entangled state (1291) goes into a mixed entangled state with more 
or less entanglement depending on the choice of the parameter a. 



^= A|2)(2| + ^|3)(3| + ^|4)(4| 
.oo^?^|l)(lK?^|2)(2| 



(29) 



5(3 + a2 
Cig) = l, Cig, 



5(3 + 
3(1 + 3a2) 
5(3 + a2) 



The concurrence of the asymptotic state C(f?oo) can be larger or smaller than C{g) 
depending on the value of the parameter a. 
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< |, i.e. the asymptotic state has less 



If, for instance, we take a = 0.5, then C(f)( 
entanglement than the initial state, and from the plot of the entanglement rate vs. the 
entropy rate (Figure 11) we can see that the conjecture (12U]) is always violated. 




Figure 10: Case 5: a = 0.5; Left: S{gt\\gt)', Right: E[gt 




0.2 0. 



line 



Figure 11: Case 5: a = 0.5; <7[gt] dashed line, |cT£;[f)i]| continuous li 

If, instead, we take for example a = 0.8, then the initial entanglement first diminishes 
and then increases again, leading to an asymptotic state with more entanglement than the 
initial one, as can be seen from the plot of the entropy of entanglement as a function of 
time (Figure 12). From the corresponding plot of the entanglement rate vs. the entropy rate 
(Figure 13) we can see that in this case the conjecture fl20l) is violated after some time. 



Figure 12: Case 5: a = 0.8; Left: S{gt\\goo)', Right: E[gt 
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0.2 0.4 0.6 0.8 1.0 



Figure 13: Case 5: a = 0.8; <j[gt] dashed line, continuous line 

Remark 4 In the last plot of the entanglement rate vs. the entropy rate, the cusp is 
due to the change of sign in the entropy of entanglement 5'(^t||^oo) and to the fact that in 
the conjecture f l20|) the absolute value of the entanglement rate cr£;[^t] is considered. On 
the other hand, all the other plots of the entanglement rate present a continuous behavior 
which reflects the fact that the entropy of entanglement does not change sign, i.e. it either 
increases or decreases monotonically. Finally, the plots of the relative entropy show its 
monotonic behavior under the action of completely positive trace-preserving maps; and to 
this corresponds a monotonic decreasing behavior for the entropy rate. 

5 Conclusions 

The time-derivative of the quantum relative entropy serves as a measure of how fast an open 
quantum system tends to equilibrium dissipating free energy under a quantum dynamical 
semigroup of completely positive maps generated by a Lindblad-type master equation. On 
the other hand, via a variational formulation, the relative entropy may be used as a pseudo- 
distance of an entangled bipartite state from the convex subset of separable states (relative 
entropy of entanglement); therefore, its time derivative can be interpreted as the speed with 
which a time-evolving state moves toward, or away from, becoming separable. 

Based on the expectation that the entanglement content of dissipatively-driven bipartite 
systems disappears asymptotically due to decoherence effects, in [TO] a conjecture was put 
forward, namely that the entropy rate, measured by (minus) the time-derivative of the 
relative entropy of a dissipatively evolving state and its asymptotic state, should always be 
larger than the absolute value of the time-derivative of the relative entropy of entanglement. 

However, beside being a source of decoherence, an environment can in some cases build 
quantum correlations that can even persist asymptotically: in this paper we have studied 
the fate of the above conjecture in the case of a Lindblad-type master equation that presents 
a rich manifold of asymptotic states that may be more or less entangled with respect to the 
initial states they emerge from. The entropy and entanglement rates have been explicitly 
calculated and numerically plotted for a class of initial states. It turns out that, when the 
asymptotic state is entangled, the conjecture is violated either at all times or after a finite 
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time; instead, the conjecture is confirmed in all cases when the asymptotic state is separable. 
The conjecture put forward in [TQ] should be thus reformulated as follows: 



<yE[Qt\ < cr[pt] ypt with ^oo separable . 



(30) 



While the asymptotic predominance of the entropy rate over the entanglement rate in the 
latter case has already been explained in [10] based on the fact that there are no entangled 
states in a suitable neighborhood of the separable asymptotic state, the truly remarkable 
fact about fl5U]) is its validity at alH > in all the cases that have been checked. 



6 Appendix A 

In order to explicitly solve the master equation ( 13T]) . one first writes the matrix A in diagonal 
form: 

/I + a 0\ / I/V2 I/V2 0^ 

A = U i l-a o \ W , U=\-i/^ V^/2 

Vool/ \o 1, 

and recasts the dissipative term in ([3]) in the form 

= 2(1 + a) - ft} + 2(1 - a) (s+ftS_ - ft}) 

+ S3ftS3-^{S3S3,ft} , (31) 

with S± := |(Ei ± zEa). Since = ^ E3, it follows that e^^^E^e"^^* = e^^"*E±. Thus, 
setting ft := e*-'^*fte~*'^*, in the interaction picture the master equation ( l3Ti) becomes 



dm = e*^*D[e-*^*fte* 

= 2(l + a)(E_ftE+-^{E+E_,ft}) + 2(1 - «) (E+ftE_ - 1{E_E+, ft}) 

+ E3ftE3-i{E3E3,ft} . (32) 

In order to solve it, it proves convenient to represent pt = Ylt j=i Pij{^)\'') {i\ with respect to 
the orthonormal basis (|6]). Indeed, using that 






v^|3) 
v/2|l) 




ll) 


= v^|3) 


fSall) = 


2|1) 


|2) 




S3I2) = 


-2|2) 


|3) 


= v^|2) ' ' 


S3I3) = 





|4) 





IS3|4) = 






(33) 
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one derives from (132|) the following equations : 

Pii = -4(l + a)pii+4(l-a)p33 , pi2 



Pi3 = -2(4 + a)pi3 + 4(l-a)p32 

P22 = -4(1 - a)p22 + 4(1 + a)p33 

P33 = 4(1 + a)pii + 4(1 - a)p22 - 8p33 

P34 = -4p34 



Pl4 

P23 
P24 

P44 



-12pi2 

-2(2 + a)pi4 
-2(4-«)p23 + 4(l + a)p3i 
-2(2 - a)p24 
, 



plus the complex conjugated equations for pj^, z 7^ j; whence 



Pi2e 



-12i 



-2(2-a)t 



Pl4(^) 
P34(t) 



— 4/ 

P34e 



(34) 



Pi2(t) = 

P24(t) = P24e" 
P44(t) = P44 • 

Of the remaining equations, two of them couple the off-diagonal terms pi3 and P32, yielding 



Pl3(t) 
P32(t) 



Pi3F,(t) + ^ii^#£^i^F_(t) 



P32 F+{t) + 



V4 - 3a2 
2(1 + a)pi3 + ap32 



a" 



F-{t) 



(35) 
(36) 



while the other three solutions couple the diagonal entries: 



PiiW 



3 + a2 



R+ Vl 



;i + a)2pil - 2(1 - «)P22 + (1 + a)V33 



;i + a)(3 + a2 



i?_(t) 



+ 



2(1 + a)pii - (1 - a)2(p22 + P33) 



3 + a2 



(37) 



P22(t) 



P33(i) 



:i + «) 



2(1 + a)pii 



a)^(P22 + P33) 



a) (3 + 



E-{t) 



[I + a)2pii - 2(1 + a)p22 + (1 + «)'P33 



3 + a2 
R + Vl-a2 



(38) 



+ 



3 + a2 

2(1 + a2)p33 - (1 - a'^){pii + P22) 



;i + afpii + (1 - afp22 - 2(1 - a2)p. 



'33 



(3 + a2)(i_a2) 



3 + a2 

where i? = pn + P22 + P33 = Piiif) + P22(^) + P33(^) is a constant of the motion and 

f cosh 2tV4 - 3a2 



-8t 



/cosh4t Vl-a2 
lsinh4tv/r-^ ' ^±W=e 



a 



\sinh 2tV4 - 3^2 



(39) 



(40) 



are quantities which decay asymptotically with t — t- +00. The remaining entries Pij{t) follow 
from complex conjugation. By returning to the Schrodinger representation, using f l5^ the 
explicit solution of (ET]) reads 



(41) 
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7 Appendix B 

In order to prove the Proposition in Section 3, let us consider the spectral decompositions 
St = (see dH])) and g^ep = Y.t=i Sj\sj){sj\. We have: 



Tr Qt log g, 



ep 



< 



^ ri{t){i\ log Qsepli) = ^ri{t)^\{i\sj)\'^ \ogSj 

i=l 1=1 j = l 

4 4 

^ri{t)\og(j^Sj\{i\sj) 

i=l j=l 

4 

'^ri{t)log{i\gsep\i) = TilgtlogUlg, 



isep\ j ) 



where the inequality follows from the convexity of log x, 

log ^ Ai Xi > ^ Aj log , Ai > , ^ Aj = 1 , Xi>0 



(42) 



and Yl^=i K^l'^i)!^ = 1- Also, we have introduced the completely positive map 



g h-> n[^] := ^ \i){i\g\i){i 

i=l 



(43) 



on the 2-qubit density matrices S{C^) that diagonalizes its argument with respect to the 
orthonormal basis ([6]). This map has the following property which allows one to analytically 
solve the variational problem (H2|) . 

Lemma 11 : S{C^) t— )■ S{C^) maps separable states into separable states. 

Proof: Given the density matrix of an arbitrary 2-qubit state in the standard basis 

/ £"00,00 P00,01 £'00,10 £oo,ii\ 

£"01,00 £oi,oi £01,10 £01,11 

£"10,00 £10,01 £10,10 £10,11 

\£ii,oo £11,01 £11,10 £11,11 / 

the action of the map 11 transforms it into a density matrix of the form 



/ £00, 



U[g] 



00 





eoi,oi+eio,io 







\ 



3f^e(^oi,io) 



01+£ilO,10 



3fte(f?oi,io) ^ 2 
V guAiJ 
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By partial transpostion [20], Il[g] is entangled if and only if |3fJe(f)oi,io)| > A/f?oo,oof?ii,ii- But 
then, the partially transposed g (with respect to the second qubit). 



9' 



( f?00,00 f?01,00 f?oo,io f?oi,io\ 

f?00,01 f?01,01 f?00,ll f?01,ll 

f?10,00 f?ll,00 f?10,10 f?ll,10 

\f?io,oi f?ii,oi f?io,ii Q\\.\\) 



Q\\ f?oi,io \^ Therefore, if ^) is separable, then also n[^] must be so. □ 



cannot be positive semi-definite, for |f?oi,io| > |3^6(^oi,io)| > ^f?oo,oof?ii,ii the sub-matrix 
f?ii f?oi,ic 
f?io,oi Q22 

Observe that ({42]) implies sup^^^^ Tr (^^^ log ^^ep) < sup^^^^ Tr (^^^ log n[^sep]) ; on the other 
hand, since 11 maps separable states into separable states, 

sup Tr ( gt log g^e-p ) < sup Tr ( gt log 11 [p.^p] ) < sup Tr ( gt log g.^v ) • 

Thus, the maximum in (142|) is attained on the subset 5^*^^ of separable qubit states that 
are diagonal with respect to the orthonormal basis, namely of the form ( |22l) with the sec- 
ond bound on the real parameters x, u, v in ( l23l) coming from the condition of positivity 
under partial transposition of matrices of the form ([7]), which is necessary and sufficient for 
separability. 

It thus follows from ( l2Tll and ( H2l) that for gt as in (IHl) the relative entropy of entanglement 
can be reduced to the computation of ( l24l) . In order to explicitly solve such a variational 
problem, we seek the stationary points of a function of the form 

f{x,y,u,v) := a log a; + dlogy + blogu + clogf + X{x + y + u + v — 1) 

with given a,b,c,d > such that a + b + c + d= 1, relative to variations of the parameters 
x,y,u,v over values achieving separable states of the form ( l22l) . Stationarity implies 



a = —Ax, d = —Xy, b = —Am, c = —Xv ; 

whence A = — 1 and a = x,d = y,b = u,c = v. However, this can be the required solution 
only if the state g in E[g] is separable so that E[g] = 0. Otherwise, the solution must lie on 
the border of the subset of separable states of the form (122|) . where the inequality in (123!) is 
saturated. From x + u + v + y = l and {u — vY = Axy, one gets 



^± = 2 ' = 2 

so that the function to be maximized becomes 



/(x,?/) = alogx + dlog?/ + 61og( ) +clog( j. (44) 
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Stationarity with respect to x,y leads to a system of two equations for the two unknowns 
X and y in terms of the coefficients a,b,c,d. From setting dx^yf{x,y) = and from the 
condition a + b + c + d = lit follows that y = x — a + d and that 



By inserting into it the values (!9l)- (fTT!) . this expression yields the separable state of the 
form (1221) which is closest to an evolving entangled state of the form ([S]). Though cumber- 
some, the resulting entanglement rate fll9p is amenable to numerical inspection. 
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